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Abstract 

Connection between the quasar luminosity function and galaxy mass function is inves- 
tigated in the framework of a phenomenological approach which relates AGN formation to 
galaxy mergers. Quasars are assumed to be short-lived, the luminosity of a quasar is con- 
trolled by the masses and angular momenta of the merged galaxies which have formed the 
quasar, and the amount of gas in them (the masses and momenta determine the quantity 
of mass which loses its angular momentum and can fall to the center). The proposed model 
can explain the shape and evolution of the quasar luminosity function, and allows us to 
estimate the parameters: the fraction of matter which falls into the center r] (which seems 
to be related to the quantity of gas in the galaxies) and k (an average density contrast in 
the regions where quasars form). The obtained values of k vary from ~ 4-7 at z = 0.5 to 
~ 1-2 at z = 2, r] vary from a few per cent at z = 0.5 to a few tens per cent or even values 
close to 1 at z = 2. In contrast to the cases considered earlier by the authors, the Eddington 
limit which, probably, can be exceeded in quasars plays an essential role. 



1 Introduction 

In present time, it has been established rather reliably that galaxy interaction (in particular, 
merging) correlates with the activity of galactic nuclei, at least, for powerful objects (see, e.g., 
reviews by Heckman, |1990|; Kontorovich, |1994| , and references therein). There is still some uncer- 



tainty for less powerful AGN, such as Seyfert galaxies. However, in the spirit of the unified AGN 
scheme (Antonucci, |1993|) we may suppose that here we also deal with interaction. 

Note, that such bright objects as quasars are a sort of markers: the quasar formation epoch 
which is often identified with the well-known cutoff in their distribution at z = Zcr ~ 2.5 (Schmidt 



et al., |1995| ; Shaver, p.995| ) may be also the epoch of massive galaxy formation due to mergers 
of less massive blocks (see Kats et al., |1992| ; Kontorovich et al., |199^)[] Note, that reality of the 
cutoff in the z-dependence of quasar density is confirmed by the counts of radio sources (see the 
review by Peakock, p.989| and the paper by Artyukh and Tyulbashev, |1996|) . 

Last data from the Hubble Space Telescope seem to confirm this point of view. Observations of 
galaxy formation of blocks (merging process) in the redshift range 2.6 < z < 3.9 allows Clements 



-'^The more so, as the merger process may be "explosive", see Cavaliere et al. (1991); Kontorovich et al. (1992) 



and Couch (|1996|) to conclude that, possibly, an epoch of {massive, K. K.) galaxy formation has 
been discovered. Observations of subgalactic blocks aX z = 2.39 (Pascarelle et al., |1995|) and their 
relation to galaxy formation were discussed in details for the rich group which was discovered 
in connection with a faint blue galaxies investigation program (see also other works of the same 
group: Windhorst et al., |1993| ; |1995| , and references therein). Recent source counts also allow to 



explain their rise to the past (for the standard cosmology and critical density) by evolution of the 
number of galaxies due to their mergers, assuming that this epoch corresponds to z ~ 2 (Metcalfe 



et al., |1995| ; see also Bender and Davis, p,996| ) 



On the other hand, direct observations of close quasar host galaxies by the HST (Bahcall 
et al., |1994| ; |1995|) gave a remarkable confirmation of the direct connection between the activity 
and the galaxy interaction and merging. In particular, in the case of PKS 2349 quasar host galaxy, 
a LMC-size sinking sateUite galaxy was discovered^ 

So, there are strong reasons to continue investigation of the relation between the activity and 
galaxy mergers (Kontorovich, Krivitsky, |1995| ) and perform a more detailed comparison between 
the observations and the phenomenological scheme which was proposed earlier and improved 
below. We assume that pairwaise mergers of galaxies, taking place due to their gravitational 
interaction and tidal forces, are the factor which triggers activity. In this approach, quasar lu- 
minosity function (LF) and galaxy mass function (MF) turn out to be related. Below we shall 
assume that galaxy MF is known and, thus, shall not make any assumptions about the mechanism 
of its formation. The activity, in turn, is controlled by mergers. 

In this work we shall found the quasar LF and the values of parameters (the density contrast 
and the parameter which determines probably the amount of gas in galaxies) for which it agrees 
with observational data in the redshift range 0.5 < z < 2; we shall analyze also the connection 
between galaxy MF and quasar LF asymptotical behavior. 



2 Discussion of the model 

The probable cause of the correlation between the galaxy interaction (in particular, merging) and 
their nuclear activity is that the interaction leads to redistribution of the angular momentum and, 
therefore, some part of matter (probably, gas) gets into the central region and gives material for 
accretion (see, e.g., Hernquist, Barnes, |1994| ). 



In the proposed earlier model (Kats, Kontorovich, 1991 ; Kontorovich, Krivitsky, 1995|) which 



describes appearing of activity due to mergers, the falling of matter to the center was assumed 
to be related to compensation of a part of the angular momentum at the merger. Below this 
assumption will be considered as a special case. According to this approach, the most important 
parameters of the problem are galaxy masses angular momenta S, and the amount of gas. 

We shall assume that the luminosity of an active galaxy formed as a result of merging between 
two galaxies is controlled by their masses 9Jti_2 and momenta Si^2; as well as the collision orbital 
momentum J: L = L(JMi, DJI2, Si, S2, J) (taking into account the amount of gas will be discussed 
below). To compute this function, a detailed theory is needed, which would deal with a very 
complicated multi-step process due to which a part of matter lose its momentum and gets into 
the center after the merger. However, we shall restrict ourselves by a simplified phenomenological 



^In some sense, it is rather the absolutely unperturbed eUiptical host galaxies, discovered by the same group, 
to be a puzzle. 
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approach. The number of active galaxies formed per unit time and unit luminosity range is, 
obviously, expressed as 

/(L) = y /(OJli, Si)/(!m2, S2)f/(0Jti, 0Jt2)F( J) X 

6{L - L{mi,m2, Si, S2, J)) dTli d^Si d^S2 d^J, 

Here /(9Jl, S) is the galaxy mass and angular momentum distribution function; f/(9Jli, 9JI2) is 
a characteristic of the probability of a merger between galaxies with masses Tli and DJI2 (in 
general, U depends not only on masses, but also on momenta, however, this dependence seems to 
be less essential and will not be taken into account in this work); -F(J) is the angular momentum 
distribution function. Thus, given /(OJl, S), we can found the rate of active objects formation, as 
a function of their luminosity I{L). 

Next, it is possible to relate I{L,t) to the active nuclei LF (j){L,t). To do it, we have to 
make some assumptions about the evolution (i.e, in our case, the light curves) of active nuclei 
forming by mergers. Thus, if we assume that the light curve is step-shaped, with the average 

may, in general, depend on L). Another possible case: if we assume that the luminosity of an 
active galaxy decreases exponentially, with the e-fold time tact, then (f){L,t) is described by the 

equation ^ — ^ (f^) ~ ^ analogue of the continuity equation in the luminosity space; 

cf. Cavaliere et al., 

formation rate (/ <^ 



985| ). In the latter case, large lifetime (tact ^ 10^ years) and small AGN 



d (I4\ 



) corresponds to the "luminosity evolution" : changing (p reflects, 

mainly, reducing luminosity of the existing objects. The case of small lifetime (tact ^ 10*^ years) 
for both equations corresponds to existence of many AGN generations, which change each other 
in the course of the Universe evolution ("number evolution"). In this paper we shall consider the 
case of number evolution (see more detailed discussion below, in section H). Assuming tact much 
less than a characteristic time of I{L,t) changing, we have for the former equation^ 

cj,iL,t)^t^JiL,ty, (2) 

and, for the latter equation 

t 

^(L,t)^^ I{L',t)dL'. (3) 

Qualitatively these two expressions are very similar: if I{L) has a power law region and an 
exponential decrease region, then the shape of (t){L) is approximately the same as the one of I{L) 
(we shall not consider the case when tact depends^ on L, though it can be easily done if there 
appear some observational data about a dependence of tact on L or some other galaxy parameters). 
So, in fact, we shall investigate the source I{L) in the equation for the LF. 

In our previous works we considered an expression for /(OK, S), which corresponds to the 
"anisotropic" momentum distribution. This distribution appeared in Kats, Kontorovich ( |1990| ) 



■^The solutions of both equations can be easily written in an explicit form for arbitrary iact, but below only 
the case of small tact will be of interest for us. It is possible also to write an explicit solution for arbitrary (not 
necessarily exponential) light curve. 

^In particular, the Eddington time does not depend of galaxy parameters. 
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as a solution of the generalized Smoluchowski kinetic equation (which describes galaxy mergers) 
without allowance for the orbital momentum, if there is some initial anisotropy. In this case, 
the initial anisotropy is amplified in the course of time, the momentum distribution tends to a 
^-function, and the momentum of a galaxy is proportional to its mass: 

/(OJt, S)^ $(971)5 (s-^), (4) 

where $(3Jt) is the MF. Distribution (^) is useful from the methodical point of view, because the 
asymptotical behavior of I{L) can be computed analytically for it. However, from the astrophysical 
point of view, isotropic distribution 

f{m,s) = ^m) ^ :^exp f-^sys^) , (5) 

which will be considered in this work, is more interesting. 

The mass dependence of an average mass to luminosity ratio for normal galaxies is rather 
weak. We shall neglect this dependence^ and take ^ ~ 10. Then the MF just coincides with the 
LF (except for the normalization). Below, in section ^, for computing I{L), we shall take the MF 
in Schechter's form 

$(mt) = $o9Jt-"e-™/^. (6) 

The index a ~ 1 for field galaxies and a > 1.25 for clusters (see, e.g., Binggeli et al., |1988| ; Loveday 
et al., 1992| ). As the integral which expresses the total number of galaxies for (|^) diverges, we 



shall assume $ = at small masses Wl < DJIq. Observational data obtained in recent years testify 
to possible steepening of $ at small masses. This steepening will be taken into account in section 

i 

We will use a rather common merger criterion: (i) minimal distance between the colliding 
galaxies is less that the sum of their radii [Ri + R2); (ii) the relative velocity at infinity v is 

R1+R2 



less than Vg = y Then the merger cross-section (taking into account gravitational 

focusing) is cr = vr 
expression for U : 



focusing) is cr = 7r(_Ri + R2)^ {I + v"^ / v"^) for v < Vg and a = for f > Vg. It results in the following 



uim„ m,) = av^l ^^^^^^^ ^ ^ (7) 



(cf. Kats, Kontorovich, |199CI| ; Krivitsky, Kontorovich, |1997| ; Cavaliere, Menci, |1997| ). Here the 



bar means an average over velocities; galaxy radius R relates to the mass as i? = C9K^ (/3 
1/3 corresponds to constant density, (3 = 1/2 to the observational Faber — Jackson and Tully — 
Fisher laws); C2 = (9/2) (371)1/2^2/^3^^, C1+/3 = 2{3TTy/^CG/vr^s, ~ {Cv^^JCy/^^-f^l For the 
function [/(9Jli, 93T2), it is convenient to introduce its homogeneity power u and exponents mi,2 
which describe its asymptotics for very different masses: 

U oc fmi^OJls', Tti < TI2, ui+U2 = u. (8) 

Obviously, for (0) mi = 0, = n = 2 if 971 < Mb and 1 + /3 if 971 > 9Jlb. Note that it is the 
parameters Ui^2 (that is the asymptotical behavior of U) to determine the asymptotical behavior 
ofJ(L). 

^Evolution of this ratio which reflects evolution of star formation determined, in particular, by mergers may be 



very important, cf. Madau (1997) 
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3 Asymptotics and the relation between indices 



Given the asymptotics of L(dJli, 9JI2, Si, S2, J), U(^Mi, DJI2) and fiVJl, S), it is possible to find the 
asymptotical behavior of I{L). In particular, the model predicts that I{L) has a power-law region, 
the slope of which depends on the slope of the galaxy MF power-law region. 

First we shall consider the simplest variant: "anisotropic" momentum distribution without 
taking into account the orbital momentum (J = 0). In this case momenta can be expressed 
in terms of masses (S oc DJl) and, so, L = L(9Jti, 07l2)- We shall assume that L(0Jti,0Jl2) has 
power-law asymptotical behavior at OJli <C DJI2: 

L(mti, 9JI2) oc 97^^9J^^^ mi < 9Ji2, (9) 

i.e., L(JMi, DJI2) can be expressed as an asymptotical power series with respect to both arguments. 
The right-hand part of is a homogeneous function of power A = Ai + A2. Sewing together 
the asymptotics for Tli ^ and VJHi ^ WI2, we shall assume below that L(9Jti,9Jl2) is a 
homogeneous function of power A in the whole range of DJli, 9JI2. Expression (|]) is analogous 
to for ?7(0Jli, 0Jl2)- Knowing Ai_2 is enough to find the slope of the power-law intermediate 
asymptotics of (0). 
Expression 

L = BAm, B = er]c\^, Am = rui + m2 - m, m = S/VGMR, (10) 

which was considered in Kats, Kontorovich ( |1991|) ; Kontorovich, Krivitsky ( |1995|) is a particular 
case of (^), corresponding to A2 = 0, Ai = A > 0. Here Am is the mass which has lost its 
momentum due to momentum compensation at the merger; e the accretion effectivity; 77 the 
fraction of Am, which gets to the central black hole; t^c the accretion time (we shall assume 



^ac = ^act); c is the light speed. Note that, though expression (|TOD was based on an oversimplified 
scheme of the origin of activity, the assumption A2 = 0, Ai > is much more general and seems 
to be rather plausible even without any connection with model (^^)- Its physical sense is that 
when a massive galaxy merges with a low-mass one, the luminosity is determined mainly by the 
latter mass. Thus, results obtained from ( pUj ) are more general than the model (pUj). In this case, 
the slope of the power-law region is determined by equations (15) and (16) from the cited above 
work by Kontorovich and Krivitsky. The opposite case, Ai = 0, A2 > (i.e., the luminosity is 
determined, mainly, by the more massive galaxy) was considered (equation (18) in the same work) 
in connection with the situation when the luminosity equals to the Eddington one L = L-^dd, and 
LEdd oc OJIh oc m^ (here Tin stands for the mass of the black hole in the galaxy center, OJl is the 
galaxy mass, X = h). The combined case L = min(i?Am, Lsdd) which also was considered there 
corresponds formally to a function L(9Jli,9Jl2) which is described by two different expressions of 
the form (P), with different Ai^2, in different regions. Asymptotics of I{L) for Ai > 0, A2 > can 
be calculated similarly to how it was done in Kontorovich, Krivitsky ( |1995| ) for Ai = or A2 = 0. 
Here we give the result, without the derivation, for completeness: 

( 1 - {u + 2 -2a)/X, k<0, 
7= I l-iu + 2-2a)/\-k/X2, A;>0, L<L(9Jto,/i), 

[l-{u + 2-2a)/X + k/Xi, A; > 0, L > L(9Jto, /x), ^^^) 

where k = —X2{u — 2q; + 2 — A)/A — a + M2 — A2 + 1. 
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For A; > 0, the plot of I{L) has a break; a more flat region change to a more steep one at 
L~L(9Jlo,/i). 

Now we shall consider I{L) for the isotropic distribution (|^) and with the orbital momentum 
taken into account. Unlike (^), for (^) with J the asymptotics cannot be determined analytically. 
The reason is that the dimension of integral (|I]) for (^) is much higher then for due to 
the 5-function in (1). In the same time, there are some heuristic arguments which lead to the 
supposition that the results mentioned above will not essentially change. Indeed, for (^) L was 
completely determined by the masses, and the 5-function in the integral cut a one-dimensional 
integration path in the (3Jli,!H2) plane (fig. 2 in Kontorovich, Krivitsky, |1995| ). In the case of (0) 
L depends not only on masses, but, averaging over momenta, we can introduce OJI2). Due 

to the scattering of the momenta, not only a one-dimensional line, but also a whole band close 
to the line will make a contribution to the integral over Tli, TI2 (after averaging over Si, S2, 
J). However, if L(3Jti,!H2) can be still described by an expression of the form (^ then one may 
expect that the asymptotics of the integral will not change. 

To verify this supposition, we carried out numerical Monte Carlo simulation. Momentum 
distribution was taken in the form (^, mass distribution $(9Jt) was assumed to be a Schechter 
function (|^) for 971 > TIq and $(9Jl) = for 971 < TIq. A root mean square momentum was 
assumed proportional to the mass in the power (3 + /5)/2 (such a dependence is formed by mergers 
in the case U oc {Til + Tl2){Mf + m^) and is close to the really observed one, see discussion in 
Kontorovich et al. (|1995|) ; Krivitsky and Kontorovich ( |1997|) ). The luminosity was calculated as 
L = min(i?Am, LEdd), where m = S/ \J GDJIR. Unlike the previous section, the orbital angular 
momentum was taken into account too. 

In this case we may expect that an effective value of A will be 1, Ai = 1, A2 = (as we 
mentioned earlier). Indeed, for DJli <^ DJI2 average proper momenta are oc 9Jl?+'^, oc 971^+^ 



oc 



the orbital one0 is oc {TliVgiRi + R2)f oc Ti\Tll , after the merger = Sf + + J 
(1 + O (il)). ^ = - ra., (1 - O (mx)). A.n ^ + oc 

We used the following simulation algorithm (a simplified description): 

1. Two random numbers, 97li,2, distributed according to the given MF $(971), were simulated. 

2. Two random vectors. Si, 2, with distribution (|^), were simulated. 

3. Galaxies 1, 2 merged with the probability proportional to f/(97li, 97I2). 

4. According to the merger cross-section assumed in our work (see page ^), the impact pa- 
rameter and the relative velocity were simulated, then the merger orbital momentum J was 
computed. 

5. The black hole mass 971h was simulated (variant 1: 971h = C^Tt''; variant 2: 97tH is an 
independent random value with a power-law distribution). 

6. Using 97li, 97l2, 9Jt = 97li + 97l2, Si, S2, J, S = Si + S2 + J, 9JIh the luminosity of the active 
object L was calculated. 



^The orbital momentum — (97119712/(9711 + 97l2))^u^p^, the impact parameter ~ (i?i + i?2)^Wg/'y^, so, 

1^ (x.m\m)^^ . 
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Figure 1: An example of the results of numerical simulation for I{L) (which determine the quasar 
LF) in the case of the isotropic momentum distribution and the orbital momentum taken into 
account. The luminosity in the figure is given in units Li (defined in Kontorovich, Krivitsky, |1995|) , 
I{L) is normalized to 1. The values of the parameters are: /x/OJlo = 10^, P = 1/3, ^Mc/VJlo = 10^, 
the exponent of the MF is a = 1.25. The plot confirms the asymptotic expressions obtained 
analytically (taking into account the hypothesis of Aefr = 1, see the text). As OJIq <C 971c ^ 
Ui + 1 — a < 0, the formula I{L) cx ^-^+('^1+^-°')/^ should be applied. For Aeg = 1 the power-law 
region of I{L) should be L~^'^^, the figure confirms it. 

Thus, the algorithm gave a random value as an output, which was distributed according to the 
same law as the desired luminosity. Repeating the computations many times, it is possible to find 
its distribution function, i.e., I{L). 

The procedure of simulation which was actually used was a bit different from the simplified 
scheme given above. In item |^, the simulated galaxies merge with some probability p; with prob- 
ability (1 — p) they are rejected. The probability p must be proportional to f/(07li, OK2), e.g., if 
U(JXHi,VJl2) is bounded above by some OTTmax, we may chose p = f/(njli, njl2)/f/max- However, in 
our case f/(07li, 97I2) is an unbounded function. Moreover, even if we introduce a limit galaxy 
mass 97Tmax and the corresponding value f/max then p = [/(OJli, 0Jl2)/t/max will be very small for 
majority of galaxies and they will be rejected. This will cause very large computation time. To 
overcome this difficulty, we transformed the integrand in (|l|) as 

j fif2Ui2S{L - Li2)F(J) £jdTli d^Si d^S2 = 

2 j mihf2^y^f^5{L - L,2)F{3)d^J ddK,ddK2d'S,d^S2 (12) 

(A = /(ajii,Si), /2 = /(!m2,S2), 

f/12 = f/(!mi,9H2), L12 = L(9Jli,9Jl2,Si,S2, J)). 

In item | of the algorithm we took f/(9Jli, 9JI2) /(9Jl^ + OJt^) instead of f/(mti, 9JI2), and in item I 
we took 9715*$ 1 and $2 instead of $1 and $2. 

The results of simulations confirmed the above supposition: the slope of the power-law region 
of /(L) coincides with the predicted value, and Aefr = 1 = Ai, A2 = 0. As an example, fig. |I] shows 
the distribution function for a particular case (see the parameters in the caption). 
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Thus, in the most interesting case, when A2 = 0, Xes = 1, Mi = 0, M2 = 1 + /? or 2, the slope 
of the power-law region of /(L) (and, therefore, the active objects LF too) just coincides with 
the slope of the galaxy MF a, according to equation (16) in Kontorovich, Krivitsky ( |1995| ). This 
result agrees well with observational data: according to Binggeli et al. ( p.988|) ; Boyle et al. ( |1988|) , 
both a and 7 are close to 1 (somewhat more). 



4 Luminosity function of active objects 

In the previous sections we described rather a general approach which relates AGN formation to 
galaxy mergers. Below we shall consider a concrete application of this approach. The purpose of 
this section is to obtain the LF from (0), (|^) and find the parameters for which it agrees with the 
observed one. 

We shall take the following input data. 

1. Galaxy mass function. The bright end of LF (and, so, MF) of normal galaxies is described 
well by Schechter's formula (|^). According to the data obtained in last years, the LF steepens 



at its faint end, the slope a reaches ~ 2 (see, e.g., de Propris et al., |1995| ; Kashikawa et al., 
1995t Loveday, [19971) . So, we shall take MF at z = clS 8b Schechter function with an 
additional break: 

$(9Jt) = $0(1 + (9Jl/9Jlbr)"'~°')9?l-"2e-^/'^. (13) 

We shall assume ai = 2, ^2 = 1-25, /i corresponds to the magnitude Mb = —21 (the mass to 
luminosity ratio being ^ ~ 10), $0 = 5 ■ 10"^/!"^"^ Mpc~^, and the break OJlbr corresponds 
to Mb = — 16. Possible change of MF with z will be discussed below. 

The momentum distribution will be taken in the form (^), with the root mean square 
momentum 



1/2 



(Krivitsky, Kontorovich, 1997 ). 



const = 0.1 (14) 



3. Merger probability. We shall use model (0), assuming C = ^ = (2.ioii'^a^)/j ; = 1/2, 
and Vrms ~ 100(1 + 2;)"^/^ km/s, which corresponds to OJlb ~ 10^°(1 + z)~^/^ TIq (such a 
dependence v{t) takes place in the linear gravitational instability theory for i7 = 1). 

4. Active galaxy luminosity. We shall assume L = min(LEdd, £^''?c^iac^Am), e ~ 0.1; the (bolo- 
metric) luminosity is calculated according to L/Lq = b ■ iq^-*(^bq-Mb) ^ ^j-^g bolometric 
correction factor being 6 ~ 10 (cf. Sanders et al., p.989| ; Weedman, |1986| ). 



Black hole mass. Correlation of the masses of black holes and host galaxies was discussed 
by Kormendy, Richstone ( |1995|) p] (for spirals, the bulge masses were taken instead of the 



galaxy masses). It was found, that, in average, DJlu oc VJl. These results were confirmed by 
recent HST data (Press release No. STScI-PRC97-01). However, there is large scattering of 



''Independently, this correlation was considered in Kontorovich, Krivitsky (1995) 
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the ratio dJiu/Tl around its average value. Thus, Tin ~ 2 ■ 10^ TIq in our Galaxy, whereas 
OJTh ~ 3-10^ DJIq in M87. Below we shall use a more complicated model than in the previous 
section: we shall assume 9JIh = where C is a random value the decimal logarithm of 
which is distributed uniformly in the range —3 ± 1. 

It is well known that the observed quasar LF essentially depends on the redshift: in the past 
quasars were much brighter than now (Boyle et al., |1988|) . There are two points of view on this 
fact in the literature. The simplest interpretation is that active nuclei are comparatively long- 
lived objects (with lifetimes of billions years), but the luminosity of each quasar decreases in the 
course of time ("luminosity evolution"). However, this hypothesis encounters some difficulties. In 
particular, if we assume that the luminosity of high-z bright quasars cannot exceed the Eddington 
limit then the contemporary active galaxies, as their "descendants" , must have very massive black 
holes (> 10^ 9^0), which seems to be ruled out by observational data (e.g., Schmidt, |1988| ). The 



other point of view is that quasars are a comparatively short (< 10* years) evolution stage of 
the majority of galaxies. Thus, Haehnelt and Rees ( |1993| ) assume that an active nuclear grows 
in almost every galaxy just after its formation, the initial luminosity equals the Eddington one 
and so the LF reflects the MF of black holes. The physical reason of decreasing the luminosity, 
proposed by these authors, is that more massive galaxies form later, their density is lower, the 
central potential well is less deep, and the black hole forming there is, in average, less massive. 
So, negative correlation between the mass of the newly formed galaxy and the initial mass of 
the black hole in its center is assumed.0. Small and Blandford (|1992|) assume that the observed 
break in the quasar LF is associated with the transition between the two modes of accretion: the 
continuous one (L = LEdd) and the intermittent one (L < L^dd and is controlled by the amount 
of "fuel"). So, decreasing of the break luminosity is associated with decreasing of the average 
amount of "fuel". Note that both Haehnelt and Rees ( |1993| ) and Small and Blandford (|1992|) 



assume that the lifetime of an individual quasar is much less than a characteristic evolution time. 

How can the evolution of the quasar LF be described in the merger model? In this model 
the characteristic luminosity corresponding to the break is related to the mass fi in (|^), (0). 
Since less massive galaxies form earlier, /i cannot decrease with time. So, cosmological evolution 
of $(971) cannot be the cause of the decreasing of the quasar luminosity. One of the possible 
explanation (which we shall assume in this work) is cosmological evolution of rj (fraction of mass 
which actually gets into the center). Here quasar lifetime is assumed much less than the age of 
the Universe, so we may use (0), (^. We shall take tact ~ 10* years. Cosmological decreasing of 
rj may be caused, for example, by decreasing of the amount of gas in galaxies. Indeed, gas and 
stars behave in different ways at merging; the matter which gets into the center is, mainly, gas 
(Hernquist, Barnes, |1994| ). 



To obtain the quasar LF for such a high z as 2, we must take into consideration cosmological 
evolution of $(9Jl). Possible reasons of such evolution are galaxy mergers, birth of new galaxies 
of gas, etc. Reliable observational data on the normal galaxy LF are available only for moderate 
redshifts [z < 0.5), and the main contribution to the change in this LF seems to be given by the 
change of the star formation rate rather than alteration of the MF (e.g.. Small et al., 1997|) . Thus, 



we have to take the time evolution of $(971) from model theoretical calculations. The existing 
theories for galaxy formation cannot yet give a detailed and reliable description of this process. 
However, they make some qualitative predictions. The evolution of the LF was considered, e.g., by 

^Siemiginowska and Elvis ( |1997D also suppose that active black hole masses decrease with time. 
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Kauffmann et al. ( 1994 ); Cole et al. (|1994 ). They concluded, in particular, that an average galaxy 



luminosity at 2 ~ 2 is several times less than at 2 = 0, whereas the amount of dwarf galaxies is 
several times higher. In the same time, both groups of authors notice that they cannot account for 
the observational data for both faint galaxy counts and the slope of the galaxy LF simultaneously: 
the calculated LF has much higher a than the observed one. Gnedin ( |1996| ) obtained the value of 
a and the shape of the MF which agree very well with the observed ones (may be, except for the 
observed steepening at very small masses) but he did not compute the time evolution of $(971). 

In this paper we shall compute (/)(L, t) for two variants of $(9Jl, t\. 1. non-evolving $(OJt) in the 
form (|T3|); 2. "maximal" evolution of $(971). In the second variant we shall take a composed MF 
and assume that the shape of $(971) is described by (0), /i and $0 depend on time approximately 
as in fig. 3 by Kauffmann et al. (|1994D and fig. 19 (right bottom) by Cole et al. (|1994| ), but. 



in the same time, a\ and 02 has the observed values (see page H), i.e., the MF is more flat 
than in Kauffmann et al. ( |1994|) and Cole et al. (|1994| )P[ Namely, we assume /x oc (1 + 2)"^/^, 
$0 oc (1 + zf!'^. The break mass OJtbr is assumed to be constant, as there are no data for its 
evolution. 

As merging occurs, mainly, in higher density regions, we should take into account the non- 
homogeneity of the galaxy spatial distribution. Due to this inhomogeneity, lijj) depends on 
coordinates, and the right-hand part of (|ip will contain a factor (p/ (p))^, where (p) is the average 
galaxy density (because the MF will be p/(p) times higher). Integrating over a large vol- 

ume V , we obtain that the average density of quasarsH is k = , .^^^ times higher as compared 
to the homogeneous situation, and k can be expressed as 

_ Jp'dV _ JpdM _ JpdN 
^ {pyv {p)M (p)iV' ^ > 

where dAi = pdV is the mass in the volume dV, dN is the number of galaxies in this volume, 

will be referred to as the average density contrast. 

We will use two 2;-dependent fitting parameters in comparing the quasar LF with the observed 
one: the fraction of matter rj and the average density contrast k. 

The results are shown in figs. |^ (the AGN LF) and § (the corresponding values of the parame- 
ters for which 0(L) has best agreement with the observational data by Boyle et al. ( |1988|) for the 
quasar LF). The figures show that the model presented above is able to account for the observed 
evolution of 4>{L). In the case of non-evolving MF the fraction of matter which gets into the center 
Tj changes from ^ 0.12-0.3 for z 2 to ~ 0.025-0.043 for z ^ 0.5, whereas the average density 
contrast in the regions of quasar formation k is ~ 1.4-2.4 for z ^ 2 and ~ 4.5-7.2 for z ^ 0.5. In 
the case of "maximal" MF evolution the parameter t] for large redshifts is much higher: ^ 0.8-1.3 
for z ^ 2 (see the discussion below), ~ 0.05-0.11 for z ~ 0.5, and k is somewhat lower: ^ 0.7-0.9 
for z «i 2, pa 3-4.6 for z ~ 0.5. Taking into account the Eddington restriction gives in most cases 
an increase of rj and k; for r] is somewhat higher, whereas k, is somewhat lower, as compared 

to m. 



^We have chosen the cases of the fastest evolution of the MF, and neglect the difference between the MF and 
the LF evolution. 

-'Here we average the equation for cf) (or (^), (||)) and keep the notation (f) for the averaged LF. 
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Figure 2: Quasar luminosity function: observed (Boyle et al., |198^ ) and predicted by the model 
described in section |. The corresponding values of the parameters r] and k are shown in fig. ^. 
Symbols a, □, v, ■ stand for observational data for z = 0.3-0.7, 0.7-1.2, 1.2-1.7, 1.7-2.2 (with 
go = 0.5), sohd lines show the merger model results for z = 0.5, 0.95, 1.45, 1.95. Figures (a)-(d) 
show the results for no Eddington restriction. The left panel corresponds to @, the right one to 
d^); Figs, (a) and (b) correspond to non-evolving mass function, (c) and (d) are the same for the 
evolving mass function. Figures (e)-(h) show the same with the Eddington restriction. Figures 
(g) and (h) do not show the plots for large z, see the text. 
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Figure 3: Mass fraction t] which gets into the center and density contrast k, necessary for agreement 
between the predicted and observed luminosity functions (fig. 0). Sohd hne corresponds to (|^), 
dotted hne to @, 1 stands for the case without the ahowance for the Eddington restriction, 2 for 
the case with the restriction, (a) is for non-evolving mass function, (b) is for the evolved one. 



Alteration of tac, b results in alteration of rj and k, according to 

7] oc t.^cS'^h, K oc (16) 



Besides, rj increases with decreasing the average angular momentum in (|1^ 

In some cases (namely, the evolving MF and (^) the values given in fig. ^ fall beyond the 
physically allowed range (the fraction of matter which falls into the center cannot exceed 1, and 
the density contrast must be higher than 1). It does not mean that the model fails: according to 
([l6|), the values falls into the required range if, for example, tac is 5 ■ 10^ years instead of 10*. 

Note that the obtained rj values for z k, 0.5 (several per cent) have the same order of magnitude 
as an average gas fraction in modern galaxies, k ^ 10 corresponds to an average density contrast 
in the large-scale structure filaments, and the value z ~ 1 corresponds to the epoch of intensive 
star formation accompanying by decreasing of the amount of gas. Next, if the parameter rj is 
really related to the gas fraction in galaxies, then such high r] values for 2 ~ 2 as several tenth 
looks rather natural: a large quantity of gas has not yet turned into stars. However, for galaxies 
with such a high gas fraction the merger criterion described in section ^ should be modified, as 
well as the expression for [/, because in a high-speed collision {v ^ fg) two stellar systems will 
pass through each other, whereas two colliding gas clouds will show quite a different behavior. 



forming a dissipative discontinuity system (e.g., Chernin, |1996[ ). 

In the variant with the "maximal" MF evolution and the Eddington luminosity taken into 
account, the obtained LF at large redshifts {z > 1.5) disagrees with the data by Boyle et al. (|1988[) 



for any values of the parameters: the number of the brightest quasars {L ~ 10'^'^ erg/s) is much 
lower than the observed one (that is why the last two curves are not shown in figs. ^ and ^h). It 
is related to the influence of the Eddington restriction. Indeed, if we assume L < L^dd, then such 
quasars must have a black hole of a mass OJl ~ 10^° ^JIq. In the same time, the "maximal" MF 
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evolution assumed here implies that the masses of galaxies at 2; ~ 2 are approximately one order 
lower than the modern ones, whereas the black hole masses DJlu = (DJl. Thus, there are too few 
black holes with 971 ~ 10^° even in spite of the ( scattering assumed here. There are many 
possible reasons for this discrepancy: 1. the luminosity may be much higher than the Eddington 
one due to anisotropy of the quasar "central machine"; 2. in a model where periods of activity 
alternate with pauses, the peak luminosity may be higher than the Eddington one; 3. gravitational 
lensing may cause an increase of the apparent brightness; 4. the case for "maximal" MP evolution 
may not be realized. Allowing for the former two factors should result in an increase of the 
luminosity (therefore, lower t] required) and decrease of the normalization (higher k required), 
which gives one more explanation of the curves in fig. ^b which falls beyond the allowed region. 

Thus, the merger model can explain the observed shape and evolution of the quasar LP and 
give an estimate for the parameters t] and k. Also, in the case of evolving MP the results agree 
well with the hypothesis that the quasar luminosity may much exceed the Eddington limit. 
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